Magnon softening in a ferromagnetic monolayer: a first-principles spin dynamics study 
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We study the Fe/W(110) monolayer system through a combination of first principles calculations 
and atomistic spin dynamics simulations. We focus on the dispersion of the spin waves parallel 
to the [001] direction. Our results compare favorably with the experimental data of Prokop et al. 
[Phys. Rev. Lett. 102, 177206], and correctly capture a drastic softening of the magnon spectrum, 
with respect to bulk bcc Fe. The suggested shortcoming of the itinerant electron model, in particular 
that given by density functional theory, is refuted. We also demonstrate that finite temperature 
effects are significant, and that atomistic spin dynamics simulations represent a powerful tool with 
which to include these. 



Recent progress in experimental techniques have al- 
lowed the first observation of the magnon dispersion spec- 
trum of a single ferromagnetic monolayer 1 . The measure- 
ments, which were performed on an atomic monolayer 
(ML) of Fe on W(110) using spin-polarized electron en- 
ergy loss spectroscopy (SPEELS), revealed magnon en- 
ergies that are much smaller compared to the bulk and 
surface Fe(110) excitations. This strong magnon soften- 
ing is in contradiction with theoretical predictions based 
on an itinerant electron model at T = K 2 3 . This dis- 
crepancy raises the possibility that a ML of Fe on W(110) 
may not be a simple itinerant ferromagnet, as generally 
assumed, and indicates a shortcoming in the theoretical 
understanding of low-dimensional systems. In fact, from 
their experimental results the authors of Ref. 1 suggested 
this possibility. Furthermore, a recent theoretical study 
by Grechnev et al has suggested that effects of elec- 
tron correlations might be more important for surfaces^. 
Given the fundamental nature of spin wave excitations, 
and their role in physical processes such as fast magne- 
tization reversaP and current induced magnetic switch- 
m £p>HH] ^ has become both timely and important to ad- 
dress these questions. 

Low-dimensional magnetic structures lack inversion 
symmetry. Consequently, significant anisotropic contri- 
butions to the magnetic ordering can arise , in t he form of 
Dzyaloshinskii-Moriya interactions (DMI) 9 10 and mag- 
netic anisotropy energies. Both of these stem from rel- 
ativists spin orbit coupling, and may be energetically 
strong enough to compete with isotropic exchange inter- 
actions, leading to complex magnetic ground sta tes in 
systems such as Mn on W(OOl) 11 and Fe on WfllOpU^I 
Furthermore, the interaction between the Fe monolayer 
and the W substrate is known to be significant, and in 
the W(001)case this can even lead to antiferromagnetic 
ordering^ElJ. An additional challenge when modeling 
spin waves is posed by temperature: the experimental 
data in Ref. [T]were taken at 120 K, which corresponds to 
a value for T/T c = 0.5, where T c is the Curie tempera- 
ture. In this regime, finite temperature effects, which are 
normally excluded from density functional theory, can be 
expected to play an important role. 

We demonstrate that atomistic spin dynamics (ASD) 



simulations based on first principles theory provide a 
powerful tool for studying magnons at finite tempera- 
ture. Using this approach, we are able to account for 
most of the experimentally observed magnon softening 
in one ML Fe on W(110)H We mainly attribute the soft- 
ening to hybridization with the underlying W substrate 
combined with finite temperature effects. We show that 
the itinerant electron model of surface and thin film mag- 
netism, as given by density functional theory, does not 
have to be abandoned, provided dynamical aspects and 
finite temperature effects are considered. 

We performed atomistic smn dynamics simulations^ 
using the UppASD package^. In the simulations, the 
isotropic exchange interactions are treated by introduc- 
ing the classical Hamiltonian 

#3 

where i and j are atomic indices, is a classical atomic 
moment and is the strength of the exchange interac- 
tion. Relativistic effects are taken into account by adding 
one or both terms from the equation 

n SO = K^2( mi - e K ) 2 + J2^J ( m * x m o) ' ( 2 ) 

where K is the strength of the anisotropy field along the 
direction of e^, and is the Dzyaloshiskii-Moriya vec- 
tor 9 -^. From this Hamiltonian, the effective interaction 
field experienced by each atomic moment is calculated 
as = —-§^7- The temporal evolution of the atomic 
spins at finite temperature is modeled by Langevin dy- 
namics, through coupled stochastic differential equations 
of the Landau-Lifshitz form, 



= - 7 m,x[B,+b,(t)]- 7 -m,x{m, x [B* + !>*(*)]}. 
at rn 

In this expression, 7 is the the electron gyromagnetic ra- 
tio, and hi is a stochastic magnetic field with a Gaussian 
distribution, the magnitude of which is related to the 
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temperature and the phenomenological damping param- 
eter a, which eventually brings the system to thermal 
equilibrium. 

We have calculated the exchange parameters, J^ , from 
first principles calculations. A supercell geometry con- 
taining seven W layers, terminated with a single Fe layer 
on each side and surrounded with vacuum, was used. 
In all calculations, the in-plane lattice constant (ai at ) 
was fixed to the experimental value of bcc W (3.165 A). 
The geometrical ground state was obtained by perform- 
ing a structural relaxation of the slab using the projec- 
tor augmented wave (PAW) methocPS, as implemented 
in the Vienna ab-initio simulations package (VASP/^4 
We found an equilibrium geometry in which there was 
a significant relaxation of the Fe-W distance (-13%), in 
agreement with previous reports^HD 

The relaxed structure was used as the input for the cal- 
culation of the Jij values, using both t he Lie chtenstein- 
Katsnelson-Gubanov method fLKG Mp^5], and the 
"frozen magnon" approximation^^. In the first case, 
the exchange interaction parameters are obtained from 
small angle perturbations from the ground state, using 
the magnetic force theorem 28 . In the second case, they 
are obtained by an inverse Fourier transform of the total 
energy E(q) for a large number of spin-spirals with wave 
vector q. For a ferromagnet, this latter method also has 
the advantage of directly providing the adiabatic magnon 
spectra at therefore providing a benchmark with 
which to compare the results of the ASD simulations, and 
to verify that the spin dynamics Hamiltonian is correct. 

In order to study the effect of finite temperature on 
the exchange interactions, we performed the LKGM cal- 
culations both for a ferromagnetic configuration, corre- 
sponding to the magnetic order at K, and in Disor- 
dered Local Moments (DLM) states^^l, which give a 
better description the magnetic order at high tempera- 
tures. The LKGM calculations were performed using a 
Korringa-Kohn-Rostoker (KKR) method, while a Linear 
Muffin Tin Orbital method (LMTO)P was used for the 
frozen magnon calculations. Both KKR and LMTO cal- 
culations were performed within the atomic sphere ap- 
proximation, using the Local Spin Density Approxima- 
tion. The anisotropy constant we used in the spin dy- 
namics simulations was taken from expe rimental data, 
and corresponds to 2K e ^S = 4.6 me\^. 

The LMTO-based spin-spiral calculations suggest that, 
rather than a ferromagnetic solution, the global energy 
minimum is given by a spin-spiral with q = [0 0.1] 
A -1 , which is in agreement with recent full-potential 
plane wave calculations 33 . This spin-spiral state was also 
found when evaluating the exchange parameters obtained 
from the LKGM calculations for the ferromagnetic con- 
figuration, but with an energy difference on the same 
scale as the maximal energy resolution of the method. 
The exchange parameters calculated from DLM configu- 
rations correspond to a ferromagnetic solution. Although 
a spin-spiral energy minimum appears to contradict the 
widely held view that 1 ML of Fe on W(110) is ferromag- 
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Figure 1: (Color online) Calculated exchange interaction pa- 
rameters Jij for a ML Fe on W(110) as a function of distance 
(in lattice constant ai a t)- The J^s labelled FMA were ob- 
tained using the "frozen magnon" approximation while the 
other curves was calculated using the LKGM method for a fer- 
romagnetic solution (LKGM), and disordered local moment 
state (DLM). Also shown are the calculated exchange inter- 
action parameters for bulk bcc Fe. The inset shows the geom- 
etry of the ML and position of neighbour j relative to site 0. 
The Jij values in the figure have been scaled with the square 
of the magnetic moment to follow the standard convention of 
Ref.Hfl 

neticLLLiLS', the difference in energy between the ferromag- 
netic and spin-spiral state is less than 1 meV. Thus, the 
ferromagnetic state is stabilized by the magnetic crys- 
talline anisotropy. This is confirmed in our spin dynam- 
ics simulations, by adding the first term of Eq. ([2| to the 
Heisenberg Hamiltonian ([!]). 

The interatomic exchange parameters calculated both 
from LKGM and frozen magnon approaches are shown in 
Fig. [I] Although the nearest-neighbor interaction domi- 
nates, longer-range interactions can also be seen to be sig- 
nificant. The obtained values display the same trend as 
recent calculations up to the fourth shell of neighbors^, 
but we note that our value for next-nearest exchange in- 
teraction is larger in magnitude. Furthermore, we have 
also calculated J^s for the unrelaxed surface structure 
(not shown). In this case, the nearest-neighbour interac- 
tion strength increases by 40%, whereas the magnitude 
of the next-nearest interaction decreases by 50%, rela- 
tive to the relaxed structure. Calculating the spin-wave 
spectra from the unrelaxed set of exchange parameters 
results in a stiffer magnon spectrum compared to the re- 
laxed system, a result that underlines the importance 
of hybridization effects between Fe atoms and the W 
substrate. We emphasize that in order to correctly re- 
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Figure 2: (Color online) Left panel: intensity of the dynamical 
structure factor S z (q, uj) for a selection of q- vectors along the 
crystallographic [001] direction for a ML Fe on W(110). Up- 
per right panel: intensity plot of the "raw" *S 2 (q, uj) directly 
obtained from an ASD simulation. Lower right panel: the 
same S^(q, cj), obtained by following the peak- finding pro- 
cess described in the text. 

produce the T = K magnon spectra from the frozen 
magnon calculation with the spin dynamics simulations, 
at least 20 shells of neighbors had to be included in the 
Heisenberg Hamiltonian. With the LKGM approach, the 
convergence was even more sensitive, with up to 80 shells 
required. 

Using the first principles exchange interactions in the 
atomistic spin dynamics Hamiltonian allows us to ad- 
dress the dy namical properties of spin systems at finite 
temperature d 19 l 34 l 35 i Of particular interest to us is the 
dynamical structure factor 5' fe (q, cj), which is the quan- 
tity probed in neutron scattering experiments of bulk 
systems^, and can analogously be applied to SPEELS 
measurements. The dynamical structure factor is readily 
obtained by a Fourier transform of the space- and time- 
displaced correlation function, 

C k (r - v',t) = (m k (t)m k (0)) - (m k r (t))(m k r ,(0)), 

where the angular brackets signify an ensemble average, 
and k the cartesian component. 

We obtain the spin wave dispersions by identifying the 
peak positions of the struct ure fact or along particular 
directions in reciprocal space^LMlMJ This approach is 
based on the adiabatic approximation, and since longi- 
tudinal spin fluctuations are neglected, the interaction 
between spin waves and the Stoner continuum is not 
taken into account. However, Stoner excitations are usu- 
ally present for short wave lengths, so the stiffness of 
the magnons should be captured by our method. Given 
the stochastic nature of the simulations, the calculated 
S k (q_, oS) spectra are very diffuse. Unless a considerable 
averaging of ensembles is performed, it can prove difficult 



to find the positions of the peaks corresponding to spin- 
wave excitations. In order to simplify the identification 
of the intensity peaks, we performed a post-processing 
scheme in which the intensity for each q vector is con- 
voluted with a Gaussian, normalized to unity, and then 
used in a power function. This allows us to decrease the 
number of ensemble averages needed for each simulation 
by a few orders of magnitude. 

Fig. [2] displays an example of the calculated dynami- 
cal structure factor 6' /c (q, a;), plotted after the convolu- 
tion with a Gaussian function, for a selection of differ- 
ent q vectors along the [001] direction of the Fe/W(110) 
monolayer. For low q values, the spectral intensity of 
S k (q_, lj) is high, and decreases with increasing q. At the 
same time, the peak width can be seen to increase with 
increasing energy. This is made clear by observing the 
peaks located at q = 0.8 A -1 and q = 1.2 A -1 . These are 
positioned at similar energies, and do not exhibit notice- 
able broadening as a function of q. The broadening of 
higher energy excitations can be explained by the damp- 
ing term in the Landau-Lifshitz- Gilbert equation. 

In our proposed approach for calculating the magnon 
spectra, temperature effects enter in two ways, explicitly 
as the temperature at which the spin dynamics simula- 
tions are performed, and implicitly by the choice of which 
magnetic configuration the interatomic exchange interac- 
tion parameters are calculated from. In order to quantify 
temperature effects, we have calculated the spin wave 
stiffness of the Fe monolayer, calculated from LKGM 
simulations starting from both ferromagnetic and DLM 
states. The spin wave stiffness was evaluated by a least- 
squares fit of the magnon curves to a second order poly- 
nomial, within the range q < 0.2 A -1 . At T = K we cal- 
culated the spin-stiffness to be 160 meV A 2 when the J^s 
were obtained from a ferromagnetic configuration, and to 
be 110 meV A 2 from the DLM set of exchange interac- 
tion parameters. When introducing temperature to the 
spin dynamics simulations, we found that, for both the 
DLM and the ferromagnetic case, the spin wave stiffness 
decreased slowly up until the ordering temperature re- 
gion, whereupon the softening became more pronounced. 
Around 120 K, the temperature at which the data of 
Prokop et al 1 was measured, we find a non-negligible 
softening of the magnons of approximately 15%, relative 
to K, for each set of exchange interaction parameters. 
Thus, a correct description of the system lies somewhere 
in between the ferromagnetic and the DLM descriptions, 
which correspond to low- and high-temperature states, 
respectively. 

In order to make a direct comparison with the SPEELS 
data reported for the Fe/W(110) monolayer systemP, we 
therefore considered a partially disordered magnetic state 
with an average magnetization corresponding to 120 K. 
At this temperature, spin dynamics simulations resulted 
in an average magnetization of 85% of the saturation 
magnetization. Consequently, we performed an elec- 
tronic structure calculation for a partial DLM configu- 
ration corresponding to this state, and subsequently car- 
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Figure 3: (Color online) Comparison between magnon disper- 
sion curves along the [001] direction for a ML Fe on W(110). 
The dots are experimentally obtained data (Ref. [TJ, whereas 
the thick purple line represents our numerically obtained data. 
For comparison, the experimental spin wave spectrum of bulk 
bcc Fe (corresponding to a spin wave stiffness constant of 280 
meV A 2 ) fRef. I37p is also displayed. 

ried out an atomistic spin dynamics simulation based on 
the obtained exchange interaction parameters. The cal- 
culated magnon spectrum is displayed in Fig. |3| together 
with the experimental data of Prokop et a/. T , and the 
experimentally determined spectrum for bulk bcc Fe^3. 
The agreement with the SPEELS data is good, particu- 
larly for wave vectors around q = 1 A -1 . Furthermore, 
the softening of the magnons in the monolayer system is 
captured. The estimate for the spin wave stiffness from 
our spin dynamics simulation at T = 120 K with partial 
DLM exchange interaction parameters is 105 meV A 2 , 
which compares well with the experimental fit of 74 meV 

A 2 . 

In our calculations, relativistic effects have only been 
taken into account in the form of the magnetic anisotropy 
energy. It has recently been demonstrated that the 



Dzyaloshinskii-Moriya interaction is also important for 
the present system, most notably by giving rise to asym- 
metric magnon spectra around the F point in the Brillu- 
oin zon d 13 * 38 l Although we have reproduced this result 
by inputting physically reasonable Dzyaloshinskii-Moriya 
vectors in the second term of Eq. Q (data not shown), 
we do not find that the DMI gives rise to any signifi- 
cant change in the overall agreement between theory and 
experiment shown in Fig. [3j 

In conclusion, we have shown that combining first prin- 
ciples calculations with atomistic spin dynamics simula- 
tions provides a powerful tool for studies of magnetic ex- 
citations in low dimensional systems. Applied to the ref- 
erence system of a monolayer Fe on W(110), the method 
captures a significant magnon softening compared to bulk 
Fe, a finding which is in agreement with recent measure- 
ments. The microsopic reason why the magnon curve of 
this system is much softer than bulk bcc Fe can be found 
in three intertwined mechanisms. Firstly, the chemical 
relaxation between the Fe and W atoms influences the hy- 
bridization between Fe and W states, which modifies the 
exchange interactions. Secondly, the magnetic relaxation 
suggests that a partial DLM configuration is appropriate 
for the experimental situation at hand, and this config- 
uration also reduces the exchange interaction strenghts. 
Finally, the dynamic effects also soften the magnon curve, 
which is captured by finite temperature simulations. 
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